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(~| ' Abstract 

In a previous paper we have presented a general scheme for the implementa- 
^ ' tion of symmetric generalized measurements (POVMs) on a quantum computer. 

This scheme is based on representation theory of groups and methods to de- 
compose matrices that intertwine two representations. We extend this scheme 
in such a way that the measurement is minimally disturbing, i.e., it changes 
the state vector 1^*) of the system to \/n|^) where 11 is the positive operator 
corresponding to the measured result. 
I Using this method, we construct quantum circuits for measurements with 

Heisenberg-Weyl symmetry. A continuous generalization leads to a scheme for 
optimal simultaneous measurements of position and momentum of a Schrodinger 
particle moving in one dimension such that the outcomes satisfy AxAp > h. 

The particle to be measured collides with two probe particles, one for the 
position and the other for the momentum measurement. The position and 
momentum resolution can be tuned by the entangled joint state of the probe 
particles which is also generated by a collision with hard-core potential. The 
parameters of the POVM can then be controlled by the initial widths of the 
wave functions of the probe particles. We point out some formal similarities and 
differences to simultaneous measurements of quadrature amplitudes in quantum 
optics. 
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1 Introduction 



The question of how to implement quantum measurements is an important issue of 
quantum information theory. Even though the standard model of quantum computers 
uses only one-qubit measurements in the computational basis at the end or during 
the computation,^ other measurements are also relevant for quantum information for 
several reasons. 

In quantum computing, models have been proposed where collective measurements 
on more than one qubit are necessary.^ In Ref. 3 a quantum algorithm is described 
which uses even more general measurements than the usual von Neumann measure- 
ments, i.e., they are not described by a family of mutually orthogonal projections but 
by a so-called positive operator- valued measure (POVM). 

In non-computing applications of quantum information theory, like future nano- 
science, it may, for instance, be useful to implement approximative simultaneous mea- 
surements of observables which are actually incompatible when measured accurately. 
An important example would be the position and the momentum of a Schrodinger 
particle. 

The implementation of generalized measurements is not trivial since this is also 
true for the smaller class of von-Neumann measurements.^ Even though it is known 
that it is in principle possible to reduce every POVM measurement to a von Neumann 
measurement on an extended quantum system it is little known so far about how 
to realize the required transformation by physical processes, particularly when the 
implementation should disturb the quantum state in a minimal way. In Refs. 5 and 6 
the class of POVMs is described which can be measured using hnear optics. Some 
special POVM measurements on low-dimensional spaces are described in Refs. 7-12. 

In Ref. 13 we have described a general design principle to implement symmetric 
POVMs on a quantum register where universal quantum computation capabilities are 
available. However, these implementations are not minimally-disturbing. 

Here we describe an extension of the theory of Ref. 13 such that the symmetric 
measurements disturb the state in a minimal way and apply it to the Heisenberg-Weyl 
group in finite dimensions. Since the latter actually defines a family of groups it is 
desirable to have an implementation which is efficient in the sense that its running 
time scales polynomially in the logarithm of the dimension, i.e., with the number of 
qubits. We show that this can indeed be achieved for Heisenberg-Weyl groups with 
a power of two as dimension. 

We then adapt the implementation scheme to the continuous situation. The cor- 
responding POVM provides a good example for a measurement where the feature 
of minimal disturbance makes sense: Given the motivation to measure position and 
momentum of a particle in order to monitor its motion one would clearly try to avoid 
disturbance as far as possible. 
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The general idea of the paper is to show that the finite dimensional circuits pro- 
vide a paradigm for the continuous variable implementation. By replacing the finite 
dimensional gates with appropriate analogues, we obtain also a possible measurement 
scheme even though it is not a priori clear how the required "gates'" could be imple- 
mented physically. However, we show that a modification of the gate sequence could 
in principle be realized by three hard-core scattering processes. This kind of idealized 
scattering in not unphysical since hard-core potential can be a useful approximation 
in many real collision processes. 

We proceed as follows. In the next section we recapitulate the definition of POVMs 
and recall a general scheme for the implementation of general measurements by or- 
thogonal measurements. In Sec. 3 we define the symmetry of POVMs and present a 
method for designing measurement algorithms for symmetric POVMs. In Sec. 4 we 
consider the implementation for two special classes of POVMs to illustrate the latter. 
Explicitly, we consider POVMs on qubits with cyclic symmetry groups and POVMs 
on d-dimensional quantum systems with Heisenberg-Weyl symmetry. In Sec. 5 we 
convert the implementation to quantum systems with Hilbert spaces of infinite di- 
mension and describe a potential realization by scattering processes on an abstract 
level. In Sec. 6 we compare this scheme to a quantum optical implementation of 
simultaneous measurements for the quadrature amplitudes. 

2 Minimally-disturbing implementation by 
von-Neumann measurements 

In this section we briefly outline a general scheme^ for the minimally- disturbing im- 
plementation of a POVM. Consider a quantum system with Hilbert space C^. A 
POVM consists of n operators Uj G C^^'^ with Uj > and n,, = Id where 
denotes the identity matrix of size d x d. A definition for POVMs on infinite dimen- 
sional quantum systems and an infinite number of results can be found in Rcf. 14. 
In Sec. 5 we use this more general definition but here we start with finite POVMs 
since we consider implementation schemes on quantum computers at first. Following 
Refs. 15 and 16 we define: 

Definition 1 (Minimally-disturbing measurement) Let{Ilj) be a POVM. Then 
a measurement is called minimally-disturbing if it changes the state vector according 



to 




given that the measurement result is j. 
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The motivation for this definition is given by a theorem in Ref. 15 stating that the 
above type of measurements maximizes the average fidehty between the input and 
the output state if the input is drawn from a uniformly distributed ensemble of pure 
states. 

The following lemma^ reduces the implementation of this kind of measurements 
to von-Neumann measurements (with Liider's projection postulate) in the standard 
basis. 

Lemma 2 (Reduction of a POVM to a von-Neumann measurement) Let P 

he a POVM with the n operators Ilj G C'^^'^. Furthermore, let C" he the Hilhert space 
of an ancilla that is initialized with |0). Then, a minimally- disturhing measurement 
of P can he achieved hy a measurement in the standard hasis of the ancilla after the 
implementation of a unitary U G ([^'^^^'^^ satisfying the equation 

n-l 

f/(|0)®|M/)) = 5^|j)®V^|vI/). (1) 

j=0 

Eq. (0) states that t/ is a unitary extension of the matrix 

M:=J2 \J) ® = (v^, v^, . . . , ^/U^f G C^"^"^ (2) 
j 

which is defined by P. In the following section we consider this extension for sym- 
metric POVMs. In some cases, the following observations help to show that a unitary 
U implements a POVM in a minimally disturbing way. 

Lemma 3 (Linear assignment of Kraus operators) Let U G C'^"^'^" he a uni- 
tary operating on a hipartite system which is initialized with |<I>) where |$) G C"" 
and 1^) G C^. Let the first component he measured in the standard hasis after the joint 
system has heen suhjected to the unitary U . Then the conditional post-measurement 
state is pure. 

Let ^c/,$,j he the Kraus operator descrihing the corresponding state change 



where j is the measurement result. For each U and j , the mapping |$) t-^ ^!7,*,i is 
linear. 

The proof is straightforward since the projected state of the composed system is 
a product state and the map given by the partial trace is linear. We find: 
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Corollary 4 (Minimally disturbing Kraus operators) Let P be a POVM with 
operators Uj. Furthermore, let U , 1$), and ^(7,#j he as defined in Lemma{^ If the 
equation 




holds for all j then U gives rise to a minimally- disturbing measurement of P. In 
other words, whenever is positive for each j, it defines a minimally- disturbing 

measurement for the POVM given by 

Uj := Afj,^ j . 

3 Implementation of symmetric POVMs 

In this section we analyze how the symmetry of a POVM can be used for the imple- 
mentation scheme of LemmaEl Here, we follow the approach of Ref. 13 where we have 
obtained a general implementation scheme for POVMs without consideration of the 
disturbance of the measurement process. This implementation scheme also relies on 
the unitary extension of a matrix that is defined by the POVM operators. It turned 
out, that the symmetry of the POVM leads to a symmetry of the matrix which can 
be exploited for the extension. In this section we show that a similar construction is 
possible for the minimally-disturbing implementation of POVMs. 
To begin with, we define the symmetry of POVMs: 

Definition 5 (Symmetric POVMs) Let a : G ([^d.xd ^ unitary representation 
of a finite group G. A POVM with operators IIq, . . . , n„_i is called {a, tt) -symmetric 
if there is a permutation representation n : G ^ Sn of the indices such that 

Here, Sn denotes the symmetric group consisting of all permutations of n objects. 

As mentioned above, the symmetry of a matrix is a useful tool for the implementation 
of POVMs. Here we define the symmetry of a matrix as in Refs. 17-19. 

Definition 6 (Matrices with symmetry and intertwining spaces) Let G be a 

finite group and a : G ^ (^mxm ^ ■ Q ^ (^nxn unitary representations. 

A matrix A G C"^^" is [a, t)- symmetric if it satisfies 

a{g)A = AT{g) 

for all g E G. We also write aM = Mr for the {a, T)-symmetry. We call the set 

Int((T, T) = {Ae C™^" : aA = At} 
of all such matrices the intertwining space of a and r. 
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The structure of the intertwining space of two representations can be easily speci- 
fied if both representations are decomposed into a direct sum of irreducible represen- 
tations of the group as the following lemma^^ shows: 

Lemma 7 (Structure of intertwining space) Let cr andr he decomposed into the 
direct sums 

cr = ^ (/m, ® Hj) and r = ^ (/„^. kj) 
j j 
of different irreducible representations kj of the group G. Then 

Int(a,r) = 0(C'"^x"^®/deg(.,)) 
j 

where deg{Kj) denotes the degree of nj. For rrij = and nj = we insert njdeg{K,j) 
zero columns or mjdeg{Kj) zero rows, respectively. 

The key observation used for the extension of the matrix M from Eq. Q to a 
unitary is that the symmetry of a POVM leads to a matrix M with symmetry. This 
is summarized in the following lemma which can be proved by direct calculation. 

Lemma 8 (Symmetry of a POVM and its matrix) // the POVM with opera- 
tors Hi, . . . , n„ is (cr, it) -symmetric then the corresponding matrix M is (cr^ ^ a,a)- 
symmetric with the permutation matrix representation crj^lg) = \'n'{g)j){j\. 

The following theorem explicitly shows how the (cr^ (g> cr, cr) -symmetry of M can 
be extended to a (cr^ ® a,cr (B i?V'i?) -symmetry of U where cr' is an appropriate 
representation and B a unitary. 

Theorem 9 (Implementation of symmetric POVMs) Let M be the matrix of 
Eq. for a {a, 7T)-symmetric POVM with symmetry group G. Let A and B be trans- 
formations that decompose aT^®a and o into irreducible representations, respectively. 
Then there is a representation o' of G such that BaB^ © a' is equal to v4(c7^ © a)A^ 
up to a permutation of the irreducible components. Furthermore, there is a transfor- 
mation W G Int(y4(cr7r © a)A\ BctB'^ © cr') which is a unitary extension of AMB'^ . 
Then 

U := A^W{B®B) 

implements the POVM for every unitary B. The unitary U is (cTtt ^ o-jCT (B Bo' B'^)- 
symmetric. 

Proof: We decompose a-,, ® a and a with the unitaries A E C'^"^'^" and B E C^^*^, 
i.e., we obtain the equations 

A{a^ © a)A^ = © kj) and BaB^ = (/„^. © kj) . 
j j 
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Therefore, the equation 



holds. Following Lemma [3 the matrix := AMB"^ has the decomposition 

j 

with Aj G C"'^''"^ and c/^- := deg{Kj). From Th. 5 of Ref. 13 it follows that BaB'' can 
be extended to A^a^, ® <y)A'^ , i.e., the representations ^4(0-7^ ® cr)y4"l' and BaB"^ © a' 
with 

i 

are equal up to a permutation of the irreducible components. We choose a unitary 
extension W G lnt(A(o"7r ® cr)^^ fieri? © a') of A^. This extension can be achieved 
by appending appropriate columns to the right side of since the matrix has 
orthogonal columns. We can write W = (A^|A^) for this extension if we denote the new 
columns by N. With this matrix we obtain for an arbitrary unitary B G (□»^('i-i)x"('^-i) 
the unitary extension 

A\N\N){B ®B) = {M\A^NB) 

oiM. □ 

As shown in the following section, the unitary W of Th. IHl can be chosen to be 
sparse for some POVMs. Within the standard model of quantum computing, this 
can be used for obtaining efficient decompositions into elementary gates for the cases 
discussed in the next section. Furthermore, there are methods known to decompose 
the transformations A and B (B B into products of simpler matrices. ^^"^^ 



4 Examples 

In this section we explicitly construct quantum circuits and implementation schemes 
for the minimally-disturbing implementation of two families of symmetric POVMs. 
First, we introduce the following notations: For m G N define tOm '■= exp{—2m/m) 
and let 

m— 1 

Xm := Yl + l)niodm)(j| G C™"'" 

j=0 
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be the cyclic shift of the basis vectors of an m-dimensional space. Furthermore, define 
the diagonal phase matrix 



rra— 1 



j=0 



and the Fourier transform 

F ■-- 



r— m-l 

V m ^ — ^ 

j,k=0 



We obtain the equalities FmXmF^ = Zm and Z^Xm = uJmXmZm which we will 
use in the following without proof. 

4.1 Cyclic groups 

Simple examples for our implementation scheme are POVMs operating on a qubit with 
a cyclic symmetry. Measurements with cyclic symmetry can, for instance, provide 
an estimation of time when applied to a dynamical quantum system. The reason 
for this is that the time evolution of a quantum system with energy eigenvalues 
being rational multiples of each other is periodic and the dynamics is therefore a 
unitary representation of SO (2). This leads naturally to the finite cyclic groups after 
discretization. 

Fix n > 2. We consider the cyclic group C„ = {r : = 1) with n elements, the 
unitary matrix representation a : C„ — > C^^^ with cr(j) = for 

Rn:=(l ^MeC^^^ (3) 



and the orbit of the vector ^/T/n{l, 1)^ G with respect to this representation of 
Cn- We have the POVM operators 

^''-^n{o u;i){l l)(o ^^0=n(^^ l) 

for j G {0, . . . , n — 1}. 

Applying the methods discussed in Sec. El we obtain the following unitary for the 
minimally-disturbing implementation of the POVM with cyclic symmetry. 

Theorem 10 (Implementation of POVMs with cyclic symmetry) The POVM 
with the operators of Eq. ^ can he implemented by the unitary 

U := ® /2)X|„(/„ ® F2)K^ G C^^^^^n 
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Figure 1: Circuit for the implementation of POVMs with cychc symmetry group 

on a qubit. The circuit operates on m + 1 qubits. On the right side the upper m 

qubits are measured in the standard basis. 



where K denotes the permutation matrix which is defined by 

K\2j) = \j) and K\2j + 1) = \n + j) . 

The rather technical proof can be found in the appendix. The idea is the following. 
The transformation A must diagonalize the cyclic shift X„ (S)/2. This can be achieved 
by (g) I2. Furthermore, we need the cyclic shift to obtain the correct order 
of the irreducible representations. The transformation B is trivial since a is already 
diagonal. The remaining transformation (/„ (g) F2)K'^ is the the sparse matrix in the 
intertwining space. 

If n is a power of 2 the ancilla system can be a qubit register and the unitary of 
Th. ^1 can be implemented efficiently as the following corollary states. 

Corollary 11 (Circuits for cyclic POVM) For n = 2"^, m>l, the unitary U of 
Th. can be implemented efficiently with the circuit of Fig. 

Proof: Since Fourier transforms can be implemented with a polynomial number^' of 
elementary gates, i.e., one and two-qubit gates, -F^®/2 can be implemented efficiently. 
Furthermore, the cyclic shift can be written as Xl^ = Fl^Z\^F2n with 

Zln = ® R2:^' ® • • • ® R2n ■ 

The unitary K is only a cyclic shift of qubits. □ 
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4.2 Heisenberg-Weyl groups 

The operators Xd and Zd are discrete analogues of translations in position and mo- 
mentum space. If denotes the d possible positions of a particle on a cyclic chain, 
the eigenvectors of Zd can be interpreted as positions eigenstates and the eigenvectors 
of Xd as eigenvectors of crystal momentum. Like for continuous quantum systems 
these observables are incompatible and it can be desirable to have approximative 
simultaneous measurements such that the result can be interpreted as a point in 
2 0?- dimensional "phase space". In Sec. 5 we discuss the continuous analogue. 

The basis of simultaneous measurements of position and momentum are POVMs 
with Heisenberg-Weyl symmetry. For all (i > 2, the Heisenberg-Weyl group is given 
by G = {Xd^Zd) and has order d^. For a positive operator ^ with tr(/i) = l/ci we 
consider the POVM with the d^ operators 

for A;,j = 0,...,rf-1. (5) 

The following theorem shows how to implement this type of POVMs. 

Theorem 12 (POVMs with Heisenberg-Weyl symmetry) Given the Heisen- 
berg-Weyl POVM P with the d^ operators from Eq. with fi := \a){a\/d for some 
state vector \a). Then P can be implemented by the circuit in Fig.\^where the inputs 
of the ancillas are given by \a) ® \a) with the complex conjugated wave function a. 
For general fi the ancilla input has to be replaced with the state vector 

d-i 

|7) ■■=VdY, v^.Jj) ® \k) e , 

j,k=0 

where \fji-^^ denotes the entry of ^JJi in the jth row and kth column. 

The proof of the theorem can be found in the appendix. In the following we briefly 
sketch the main points of the proof. For the decomposition of cTtt ® ex we observe that 
the permutation vr given by the action of the Heisenberg-Weyl group on the operators 
is a translation in the finite plane (Z/ciZ)^. This translation is diagonalized by the 
inverse Fourier transform -FJ ® Id at the end of the circuit in Fig. |21 This trans- 
formation already block-diagonalizes a-^^ ® a. However, the irreducible components 
are only equivalent, but not equal, to a. We apply the controlled Z and controlled X'^ 
operations to obtain equality. Hence, these operations correspond to the matrix A. 
The matrix B is trivial since a is an irreducible representation. The unitary extension 
W used in Th. is decomposed into two components. One component is given by 
the first two gates of the circuit in Fig. |2l the other is absorbed into the preparation 
procedure for the initial state. 
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la) 
I a) 



F, - 



Figure 2: Schematic circuit for the POVM with Heisenberg-Weyl symmetry. The 
basis state \j) of the control wire causes the implementation of the jth power of 
the controlled operation. On the right side of the circuit the two upper systems are 
measured in the standard basis. 



As already stated, we can efficiently implement F^ by elementary gates on a qubit 
register. We also obtain efficient implementations of controlled X and Z gates by 
concatenations of controlled Rd gates as defined in Eq. Q. Hence, we can implement 
the POVM with initial operator /i = \a){a\/d efficiently if the same is true for the 
preparation of the states \a) and \a). 



5 Continuous measurements 

Here we want to address how to implement Heisenberg-Weyl symmetric POVMs for 
continuous quantum systems such that we have also minimal disturbance. For a 
detailed mathematical description of such POVMs we refer also to Refs. 22, 23. 

The continuous degree of freedom can either be a Schrodinger wave of a quantum 
particle moving on a line (where the Heisenberg-Weyl group formalizes translations 
in position and momentum space) or a quantum optical light mode (where the trans- 
lations shift the quadrature amplitudes). The most natural representation of the 
Hilbert space of a particle in one dimension is H := L^(]R), the space of square in- 
tegrable functions over the real line. For a light mode, it is often more appropriate 
to choose the isomorphic Hilbert space /^(Nq) of square- summable sequences. In this 
section, we will focus on Schrodinger particles since simultaneous measurements of 
quadrature amplitudes in quantum optics have already been implemented.^'* We will 
compare our implementation to the latter in the next section. 

We first describe the continuous analogues of the "gates" in Fig. |21and show that 
their concatenation leads indeed to a correct implementation. Later we will discuss 
a modification of the scheme which can be implemented by hard-core scattering pro- 
cesses. The description below refers to the Schrodinger representation where the 
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position operator X, defined on a dense subspace of H, is the multiplication operator 

Xip^x) := xip^x) . 

The momentum operator is 

Pijj{x) := —i—ip{x) , 

where we have chosen the units such that h = 1. Following Sec. 3.4 of Ref. 14 (with 
a slight modification of the sign) we introduce a family {Ug^t) of unitaries 

{Us,Mx) :=e-"^^(x-t), 

which formalize shifts in momentum and position space. These unitaries define a 
measurement by the positive operators 

where \a) G 7i is a wave function which is sufficiently localized in momentum and 
position space. The probability density for the result (s, t) is tr^pUs^t) if the system 
state is described by the density operator p. The outcome (s, t) is interpreted as 
momentum s and position t of the particle in a "coarse grained phase space". We 
can clearly generalize the POVM above by replacing |a)(a|/(27r) with any operator 
p, having trace l/(27r). 

In agreement with the discussions of finite POVMs in the preceding sections we 
want to implement the POVM in such a way that the state changes according to 

^ ^ tr(n3,tp) ' 

given that the measurement outcome is {s,t). 

Now we describe how to find a continuous analogue of the circuit in Fig. |21 The 
system Hilbert space (C^)®^ is replaced by H^^, i.e., in additional to the particle to 
be measured one uses two particles in one dimension as ancilla system. The final von 
Neumann measurement is a position measurement on both ancillas'^. 

The continuous analogues of the required gates are as follows. The discrete Fourier 
transform (whose inverse is occurring three times in Fig. |2|) is replaced with the 
continuous unitary Fourier transform 

{F^){x) = y^/" e-'^ymdy. (6) 

''One could also use the remaining two dimensions of a particle in three dimensions as ancilla 
system. 
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The controlled cyclic shift is replaced by a unitary Y describing controlled translations 
on the real line. It acts on the wave function ip of two particles according to 

{Yip){x,y) = ilj{x,y -x), (7) 

since this transformation would correspond to the transformation 

|x) ® \y) 1-^ \x) ®\x + y), 

if such position eigenstates |x) and \y) existed. Conjugating Y with the Fourier 
transform on the second tensor component makes it more apparent that it is indeed 
a unitary map since we obtain then the multiplication operator 

V := (^{I ® F)Y{I ® F^)iIj^ (x, y) = e-'^'y-ipix, y) . (8) 

Here / denotes the identity operator on Ti. The unitary in Eq. (jH)) is the straightfor- 
ward generalization of the controlled phase-shift operation that is the fourth gate in 
Fig. 121 The following theorem shows that the above described replacements provide 
in fact the desired measurement procedure: 

Theorem 13 Replace the gates in Fig. with their continuous analogues as follows: 

1. Set the inverse of the continuous unitary Fourier transform given by Eq. (0j 
instead of F\ 

2. Set the inverse of Y given in Eq. instead of the controlled-X\ gate. 

3. Set V as given by Eq. (0) instead of the controlled-Zd gate. 

Let fi be an arbitrary positive operator with tr{fi) = l/(27r) and the two ancilla 
systems be in the state 7 with 

00 

I7) :=5Zv^l«.-)® (9) 

j=0 

where \aj) is an eigenvector basis of fi such that 

271 fi = ^Xj\aj){aj\ . 
j 

Then the resulting transformation on Ti®^ implements a minimally- disturbing mea- 
surement for the POVM 

n,,t ■.= Us,t^iUlt. 

when followed by position measurements on both ancillas at the end and interpreting 
the position of the first particle in Fig. \^as t and the position of the second as s. 
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Proof: Due to Lemma 01 and its corollary it is sufficient to restrict the attention 
to rank-one operators fi := |a) (Q;|/(27r) and show that the unnormalized output state 
coincides with the desired state. The linearity argument holds also if /i is an infinite 
series since one can check that the map 

is continuous with respect to the topologies induced by the Hilbert space norm and 
the operator norm, respectively. This is seen from 

||$||^ = tT^Alj^^jAjj^^j) < \\A\j^^jAu^^j\\ = 

The whole "circuit" creates some wave function ^ G H^^. After measuring t and 
s we obtain an unnormalized conditional state vector given by the wave function 

We want to show that it satisfies 

\^t,s) = V^l^) = ^J^Us,t\c^){a\Ul^) . 
This means explicitly that 

^^^^(z) = J ^e~'''a{z -t) [ a{u - t)e''''^{u)du . (10) 



2vr 

Now we calculate the effect of the circuit starting with the joint state 

where \1/ is the wave function of the measured particle. First, we apply the controlled 
inverse translation and obtain 

a{x)a{y + . 
The inverse Fourier transform changes this state to 

The second controlled inverse shift followed by the controlled phase yields 
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After applying the inverse Fourier transform to both ancilla registers we obtain 



oo POO POO 



— oo •/ — oo J — oo 



V 87r3 

We simphfy this term into 



e'^'"e'""e'""a(u)a(w + u)e-'""^(z + w)du dv dw , 



OO POO poo 

, / / e'y'^e'^''-'^''^^a{u)a{w + u)'^{z + w)dudvdw. 
The integral over v is only non-vanishing for x — z + u = 0. Hence, we obtain 

— / e''"'"a{z ~ x)a{z + w - x)'^{z + w)dw . 

With the substitution u := z + w we get 

T 





^^e-'y'a{z-x) J e'y''a{u-x)^{u)du. 

The conditional state given that we obtain the result x = t and y = s coincides with 
Eq. (HHl). □ 

In order to realize the transformation in Th. El by a physical process we first ob- 
serve that scattering processes realize quantum gates which are close to the controlled 
phase shift in Eq. ((7j): Consider two particles interacting with hard-core potential, 
i.e., the interaction energy is zero whenever their distance is larger than some a > 
and infinite if the distance is smaller than a. In Ref. 25 we have discussed the state 
change caused by such a scattering provided that the considered time scale is small 
compared to the time scale on which the width of wave packets grows by dispersion. 
We will first explain the scattering process in momentum space since the change of 
momenta of classical particles provide a good intuition about the quantum case. The 
momentum p2 of the light particle obtains a sign change since it is reflected. Due to 
the conservation of total momentum, the heavy particle acquires an additional mo- 
mentum 2p2. The vector of momenta of both particles is therefore changed according 
to a linear transformation N given by 



Pi] ^ 1 Pi + 2P2 



Neglecting irrelevant translations in position space, the corresponding linear transfor- 
mation M in position space is already given by the requirement that the 4x4 matrix 
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Figure 3: Correspondence between gates and a scattering of two particles with ex- 
treme mass ratio. We shall call the scattering a "controlled SC gate" . 
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Figure 4: Circuit equivalent to the circuit in Fig. |2l 



transformation M (B N acting on the two positions and the two momenta has to be 
symplectic. We have therefore M = {N'^)~^ and obtain in agreement with Ref. 25 



M 



1 

2 -1 



The scattering process S acts therefore on the wave function in position space by 
multiplying the coordinate vector with M, i.e., 

(5^)(x) := V(Mx), xGR^ 

We obtain 

iS^/j){x, y) = ipix, -y + 2x) . 

In order to understand the relation to the gates in Th. we may represent this 
operation by the circuit in Fig. El The "reflection" gate R corresponds to a change 
of the wave function according to 

iRip)ix) ■.= ip{-x). 

Elementary calculations show that Fig. |21 is equivalent to the circuit in Fig. 0] 
where we have absorbed the Fourier transform on the wire in the middle by replacing 
an X- measurement with a P-measurement. 

After we have converted the desired circuit into the equivalent one in Fig. |3 that 
avoids controlled phase gates we still have the problem that it requires controUed-X 
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gates and its inverse instead of a controUed-X^ gate which implements the shift twice. 
However, we observe that we may convert these gates into each other by conjugating 
them with the unitary squeezing operator 

S2\x) := \2x) 

combined with reflections when needed.*' We will see later that we do not have to 
worry about the physical realization of 5*2 since we need this gate and its inverse only 
at the end or at the beginning of the first or the second wire. Hence, they can either 
be absorbed into the preparation procedure or into the measurement by reinterpreting 
the result. 

We will furthermore modify the entangling operation on ancilla 1 and 2, i.e., the 
first gate of the circuit in Fig. 0] for the following reason. An important feature of 
the circuits in Figs. El and 0] is that the POVM consists of rank-one operators if the 
input ancilla state is the product state \a) ® |a) and entangled inputs lead to POVM 
operators of higher rank. The preparation of these entangled states was not considered 
in Subsection 14.21 Here we also want to describe how to entangle ancilla 1 and 2 when 
POVMs of higher rank are desired. The goal is therefore to change the operation on 
ancilla 1 and 2 preceding the interaction with the system to be measured such that 
a family of product input states allow the implementation of POVMs of higher rank. 
In other word, we want to tune the achieved information and the caused disturbance 
of the measurements by plugging different product states into the circuit. 

After subsequently replacing the gates in Fig.^with scattering processes combined 
with squeezing operations and reflections and modifying the entangling operation 
between the ancillas, we found that an interesting class of POVMs can indeed be 
implemented by three scattering processes as depicted in Fig. El when the initial 
ancilla states are Gaussian wave packets. 

To understand the effect of the "circuit" in Fig.j^lwe shall compute a 3 x 3-matrix 
which describes the effect of the whole circuit on the three position coordinates. For 
doing so, we recall (see Fig. El) that the masses of the particles satisfy 

mi <^ ms < 777,2 . 

First, we implement a collision between particle 1 and 2. Here, the position of particle 
2 controls the position of particle 1. In analogy to the remarks above we describe the 
scattering and reflection by matrices that acts on the vector of position coordinates 
of the three particles. The scattering processes with the pairs (2, 1), (2, 3), and (3, 1) 

''Note that the described reduction of controUed-S'C to controUed-X is also possible in finite 
dimensions d. The definition R\x) | — x) is always possible and \x) i-^ \2x mod d) is bijectivc if d 
is odd. Then the ring TLjdTL allows division by 2. 
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Figure 5: Implementation of minimally- disturbing simultaneous measurement of po- 
sition and momentum by three scattering processes and one reflection. The masses mi 
and 777.2 of the two ancilla particles are extremely small or extremely large compared 
to the mass of the particle to be measured. 



correspond to the matrices 
^21 := I 1 I , ^23 := I 1 I , ^31 







The scatterings are followed by a reflection of the 2;-coordinate (the R gate). Taking 
into account that we have to concatenate the effect on the coordinates from the left 
to the right, the complete transformations in position coordinate space is given by 

^ := 'S'31 5*23 5*21 R, 

Let the initial state of the three particles be given by the wave function 

a{x)(3{y)i){z) . 
After subjecting the arguments to A we obtain 

a{x + 2y + 2z)l3{y)i,{2y + z). (11) 

In order to reduce a momentum measurement on the second wire to a position mea- 
surement we apply a Fourier transform to the state (jllj) and obtain 

- / a{x + 2w + 2z)(3{w)i){2w + z)e-''"ydw. 
2vr J -00 

With w := 2w + z we get 

^ / a{x + w + z)(3{{w-z)/2]i){w)e'^^'"-'^dw. 





TT 
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We define the integral kernel 




K,y{z, w) := \l ^a{x + w + z)(3{{w- z)/2 ) e'^^^-'^ . 



It defines for fixed x, y an operator K^^y on H by 

POO 



Note that the Kraus operators K^^y describe the unnormalized output state K^^ylip) 
of particle 3 given that we have measured x and y on the first and second particle, 
respectively (in straightforward analogy to the Kraus operators in Lemma El for the 
discrete setting). 

Now we show that Kj.^y can be obtained by subjecting Kq q to the Heisenberg-Weyl 
group elements by 

Kx,y = U.x/2-yKofiU[^^2_y . 

To see this, we observe that the translation by —x/2 in position space changes the 
integral kernel kofl{z, w) into fco,o(^ + x/2,w + x/2) and the additional translation in 
momentum space by —y/2 changes it into 

ko^Q^z + x,w + x)e'^^'"~'''> = k^^y{z, w) . 
This shows that the process in Fig. El implements a measurement for the POVM 

Hs,* := Us,tKQ QKo^oUs,t , 

when reinterpreting the measurement outcomes x,y on particle 1 and 2 as t = —x/2 
and s = —y, respectively. In order to obtain a minimally-disturbing implementation, 
we have to ensure that Kq q is positive (in straightforward generalization of Corollary|31 
to the continuous setting) because it can then be interpreted as ^/JI. If a and /3 are 
real and P is an even function, i.e., (3{—y) = P{y), Kq q is self-adjoint due to 



ko,o{z,w) = kofi{w,z) . 
The integral kernel of Kq q is explicitly given by 

kofi{z,w) = -^=a{w + z) p(j_^w - z)/2^ 



. 27r 

Now we assume that a and (3 are both real Gaussian wave functions with widths ai 
and (72, respectively, i.e., 

, , 1 / \ 1 ^/ X 1 f 

a[x) := . exp - — ^ and I3[y) := , exp 
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Under these conditions, fco,o defines a positive operator whenever ai > 2a2. This 
follows from the following lemma after replacing a and b with l/{2al) and l/(8cr|), 
respectively. 

Lemma 14 The operator given by the integral kernel 

k{x,y) := rfe-''(^+^)'-^(^-^)' 
with d > is for all b > a > positive. 
Proof: Rewrite the kernel as 

k{x, y) = d e-^--\-^^-)^-y)' e-^-y' . (12) 
It is known that the integral kernel 

k{x,y) ■=de-<''-y^^ 

defines for all positive c, d a positive operator^^ which we shall denote by K. Then 
the operator K given by the kernel (|T^ can be written as = DKD where D is the 
multiplication operator 

iD^lj){x) :=e-2"">(x). 

Hence, K is also positive. □ 

Since we have now described sufficient conditions for which Kq q is positive, we 
would like to better understand the POVM operator /i = K^q. As simple computa- 
tions show, it is (up to the normalization factor 27r) given by the reduced state of one 
particle in a two-particle system, if the latter is described by the wave function 

0(x, y) := v^fco.o = a{x + y)p({x - y)/2) . (13) 

It can be obtained from the state \a) ® by a linear mapping of the wave function 
arguments according to 

(:)-(v. -vO(:)-K:)' 

Such a linear operation transforms the initial Gaussian state into a entangled Gaussian 
state. Since Gaussian states are completely determined by their covariance matrix^^ 
we will compute the latter for the state in Eq. [121 
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For doing so, we must describe the linear transformation corresponding to G that 

acts on the arguments of the wave function in momentum space. According to the 
remarks at the beginning of this section, it is given by the transposed inverse: 



Px 

Py 



The covariance matrix of a two-particle state p consists of the entries 



Px 

Py 



tripXiXj) - tr(pX,)tr(pX,) with i, j = 1, . . . , 4 

where Xi,X2 denote the position operators and X^jX^ the momentum operators of 
particle 1 and 2, respectively. For the state \a) ® \P) it is given by (see Ref. 27) 



a — 



\ 



0^ 










oll2 








l/(2aD 










l/(2ai) } 



If the coordinate vector in the position wave function is subjected to some area- 
preserving linear map G and the coordinates of the momentum wave function to 
(G-^)"^, the covariance matrix transforms in the following way: 



G- 







a 







\ 








al - Aal 
al + Aal 







G 







4 I 1 









4 1 
4 I 1 



I 



as simple computations show. The covariance matrix of the reduced state of each 
particle is given by the 2x2 sub-matrices that refer to its position and momentum. 
Due to the symmetry of our state, it is for both particles given by 



1 



al + Aal 








4 , 1 
~7 ~r 



(14) 



It is known^^ that such a state is pure if and only if the determinant is 1/4. This is 
given for a\ — 2ai. One can rewrite a Gaussian state of a single mode having diagonal 
covariance matrix as a thermal state of a harmonic oscillator^ with frequency a;, mass 

'^In quantum optics, one would also need squeezing transformations to obtain a general diagonal 
Gaussian state. But here the product of frequency and mass of the oscillator provides an additional 
free parameter. 
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m and average phonon number N. It is explicitly given by 



n=0 

where \n) with n G Nq denotes the nth energy eigenstate of the oscillator. We will 
first use dimensionless position and momentum variables 

X' := ^{a + a^) = ^/^X (15) 

and 

P':=-l=(a-at) = ^P, (16) 

with creation operator and annihilation operator a. In these coordinates, the 
covariance matrix of the thermal state with average phonon number is the identity 
matrix times {N + l)/2, this follows, e.g., from Eqs. (2.16) in Ref. 27. In natural 
units, we have therefore the covariance matrix 



N+1 f ^ 



2 V ^ 

Comparing Eq. (fT^ to Eq. (fTTj) we obtain 

vmuj) = 5 — 5— 

a? + 4ai 

and 



(17) 



Hence 




£- J i 1 



For 0"! = 2(72 one obtains = 0, i.e., the ground state of the oscillator which 
corresponds to a rank-one operator fi. We rephrase the findings implied by the above 
discussion theorem: 

Theorem 15 Given three particles such that their masses satisfy 

nil ^ ^3 ^ ^2 • 
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Let the first and the second particle be in Gaussian states with real wave functions 
such that their widths satisfy ai > 2a2- Then the sequence of scattering processes de- 
picted in Fig.\^implements the Heisenberg- Weyl POVM {Ug^t Ul^)s,t in a minimally- 
disturbing way when the position of particle 1 and the momentum of particle 2 is 
measured and the result {x,y) is interpreted as t = —x/2 and s = —y. The initial 
operator fi of the POVM is given by 

1 

ZTT 

where pN,muj is the thermal equilibrium state of a harmonic oscillator with mass m 
and frequency u when the temperature is chosen such that the average phonon number 
is N. The parameters N and mu are determined by the widths ai and (J2 according 
to 

.2 4/^? + 



[muj] 



and 




1. 



We want to briefly explain qualitatively how the measured POVM is tuned by the 
parameters cxi and (T2. The ratio of both determine the purity of /z, for ai = 2a2 we 
obtain a rank-one measurement. By increasing or decreasing both we can achieve a 
better resolution in momentum space or in position space: Small values cri,a2 lead 
to good position measurements for the cost of having large errors in the momentum 
measurement. If cti 3> 2a2 both position and momentum measurements are bad and 
we obtain a measurement with small disturbance. 

For detailed discussions on the disturbance and accuracy of the measurements we 
refer also to Refs. 22,23. It is shown that the outcomes for position and momentum 
in POVMs of the above type satisfy the inequality AxAp > h in contrast to the 
Heisenberg uncertainty relation AxAp > h/2. 



6 Comparison to quantum optics implementations 

There are meanwhile several methods known to measure the quadrature amplitudes 
of a light mode simultaneously (see Ref. 24). We will consider the scheme shown 
in Fig. ini which has some nice similarities to the continuous analogue of our circuit 
in Fig. 121 In the following, we will use the dimensionless formal position and mo- 
mentum operators as given in Eqs. (|T5|l and ()16|1 which generate the momentum and 
position translations in the Heisenberg- Weyl group and define furthermore implicitly 
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a Schrodinger representation of a single mode state as wave function ip G L^(]R) in 
position state. 

The method in Fig. IHl uses also two ancilla modes. As in our proposal, the en- 
tanglement of the two modes tunes the POVM operator. One part of an entangled 
two-mode state (wire 2 and 3 in Fig. ^ interferes with the input state (wire 1) in 
a beam splitter. One of its output modes is subjected to a position measurement, 
the other to a momentum measurement. The results of the measurement determine 
furthermore displacements performed on the second component of the entangled in- 
put. The idea behind the scheme is to perform a teleportation using a non- maximally 
entangled bipartite state (a maximally entangled state does not exist anyway in con- 
tinuous variables) as resource. Then the transfer of quantum information is not per- 
fect but the measurements performed during the "bad" teleportation provide some 
information on the input state. Similar to our scheme, the more entangled the joint 
state, the less information provides the measurement and the less disturbance on the 
output state will be observed. 

One technical difference to our scheme is that the input and output are not on the 
same wire. The main difference is, however, that the interaction between input and 
entangled ancilla is given by a beam-splitter, whereas we use scattering processes. 
This is geometrically the difference between a rotation or a shear in coordinate space 
(for details see Ref. 25). Note, however, that the effect of the controlled displacements 
of Fig. ini could be mimicked by a controlled-X gate from wire 1 to 3 and a controUed- 
Z gate from wire 2 to 3 if the latter gate is conjugated by a Fourier transform on wire 
2. The reason is that it does not make a difference whether the controlled operation 
is performed before the measurement or afterwards. The scheme contains therefore 
quite similar elements as ours. 

In Ref. 24 we did not find an explicit remark saying that their implementation 
scheme is minimally-disturbing in the sense considered here since the authors use 
the term "minimally-disturbing" in a different sense. Furthermore, the attention was 
restricted to Gaussian states for both the input as well as for the ancilla states. We 
have observed that the implementation is also apphcable for non-Gaussian states and 
non-Gaussian POVM operators: 

Theorem 16 The scheme ofRef.\^can in principle be used for a minimally- disturbing 
implementation of any Heisenberg-Weyl symmetric POVM 



Us,t /i ult 



by preparing the ancilla state 




24 



input 




X measurement 



P measurement 



Figure 6: Measurement scheme of Ref. 24. The box at the beginning of mode 2 and 3 
indicate the entangled input on these two modes. The entanglement tunes the POVM. 
In the limit of infinite entanglement the output coincides with the input state and no 
information is gained. If mode 2 and 3 start in a product state, a rank-one POVM is 
implemented. The input interferes with mode 2 in a balanced beam-splitter were one 
output mode is subjected to a position measurement and the other to a momentum 
measurement. The results determine the displacements in position and momentum 
the output is subjected to. 
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where the \aj) denote the eigenvector basis for fi and Aj/(27r) the corresponding eigen- 
values. 

Proof: Due to the linearity argument in Lemma El we may prove our statement 
for the case that the two ancillas are in a product state. The initial three mode wave 
function is then given by 

<Pix,y,z) = ij{x) a{y) (3{z) , 

where ip is the wave function of the mode to be measured. The beam splitter transfers 
it to the wave function 



^((x + y)/V2^ a((a; - y)/V2^ Piz) . 

We simulate the momentum measurement by an inverse Fourier transform followed 
by a position measurement. Conditioned on the measurement result {x,y) we obtain 
therefore a one-mode wave function (having z as argument) which is given by 

j ^(^{x + w)/V2^a(^{-x + w)/V2^e''"ydw'j (3{z). 

After obtaining the measurement results x and y on wire 1 and 2, respectively, the 
conditioned displacement of position and momentum by \/2x and \/2y, respectively, 
leads to 




27r 



^(^{x + w)/V2^ a(^{~x + w)/V2j e''"ydw^ (3{z - V2x) e"'^y 




^1 J '^{w)a{-V2x + w)e''^^ydw^ (3{z-V2x)e 



This shows that the unnormalized state vector of the third particle, given that x, y 
was measured, reads 

After taking into account that the quantum optics convention for position and mo- 
mentum differs from the canonical definition of Eq. ()15|1 and (|T6|l by the factor -\/2 
(see Ref.^^), this is exactly the desired output state. 

By choosing the input |a) ® |a) we have therefore /i = |Q;)(a| = ^JJi. Similarly we 
can obtain operators with higher rank by choosing entangled input states. □ 

Note that the calculations which show that the scheme does indeed implement a 
minimally-disturbing POVM is very similar to the calculations in Sec. El which shows 
the close formal analogy of both methods. 
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7 Conclusions 



We have presented a general scheme to implement minimally-disturbing symmetric 
measurements by quantum circuits. By applying it to the Heisenberg-Weyl group, we 
obtain circuits for simultaneous measurements of position and momentum of a particle 
moving on a discrete cyclic chain. We show that an infinite dimensional generalization 
of this circuit leads to a well-defined measurement process on a Schrodinger particle 
moving on the real line using two probe particles. The "circuit" for this continu- 
ous variable quantum system can in principle be obtained by particle collisions with 
hard-core potential. The whole measurement process on the three particles shows 
some analogies but also differences to simultaneous measurements of the quadrature 
amplitudes in quantum optics using two ancilla modes. 



This work was supported by Landesstiftung Baden-Wiirttemberg gGmbH (coop- 
eration of the projects AZ 1.1322.01 and AZ 1.1422.01). 

Appendix 
Proof of Th. Uni 

The (2n x 2)-matrix M of Eq. (j2I) has the (cr,r ® o", a)-symmetry that is defined by 



as straightforward computation shows, i.e., we have {a.„ ® = ® Ri^ and 

a{j) = Ri^. To find the diagonalizing operations A and B of Th. El we observe 
that a is already decomposed into irreducible representations and B is therefore 
trivial. To decompose cTtt (S> ex we diagonalize cr,r by the Fourier transform F„. The 
eigenvalues of {FnXnF'^) ®R fi are in the order 1, uJni ■ ■ ■ i ^ i 

We 

apply the cyclic shift to group them into a sequence of pairs (07^,0;^) as in 
Lemma m Therefore, we have A = X2„(-F„ /2)- Following Th. IHlwe only have to 
find W e lnt{A{an ® (t)A\ BaB'^ © a') which is a unitary extension of N := A^MB. 
Hence, we can choose 



(X„ ® Rn)M = MR, 
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since the first two columns of tliis matrix coincide witli N. Tliis is verified by straight- 
forward computations, too. One can also easily check that W can be written as 
W = {In® F2)K^ where K is defined in Th. HUl 



Proof of Th. m 

To define M as in Eq. we have to define a correspondence between ancilla basis 
states and POVM operators. Since our ancilla system is a tensor product of two 
(i-dimensional systems, this correspondence is canonical and we obtain 

d-l 

M=J2 \J) ® 1^) ® Z'd^'dVi^X-^ Z,' e €''>^'. (18) 

j,k=0 

The symmetry {a^^ ® a)M = Ma of M is defined by 

{Ia®Xd(S)Za)M = MZd and (X^ ® ® X^) M = MX^. , 

as straightforward computations show. Following Th. IHl we decompose the represen- 
tation on the left side into a direct sum of irreducible representations. First of all, 
we diagonalize the shifts X^ in the first and second tensor components by the Fourier 
transform. We obtain 

{Id ®Zd® Zi) {Fa ®Fd®Id)M= {Fd ®Fd® Id) MZd 

and 

{Zd ®Id® Xd) {Fd ®Fd®Id)M = {Fd ®Fd® Id) MXd. 
The matrices on the left side can be written as 

(/, (^Zd®Zd) = ^ uj^^^'^'^Zd and {Zd ® ® X,) = uj\^^'Xd. 
j=o i=o 

Therefore, the representation is decomposed into a direct sum of representations that 
are equal to a up to phase factors. We now eliminate these factors. To simplify 
notation we define the block diagonal matrices 

X^od:=0Xr°'' and Z^,^ := ^ Z'/"" . 

j=0 j=0 

Using Z^XdZd = uj^^Xd and X^Z^X] = u^^Zd we obtain 

V j=0 J j=0 
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and 

\ j=o J j=0 

Using both equations we can write 

and 

X^od^L ®Id® Xd) ^divX^od = ®Id® Xd) , 
where we have no phase factors. Consequently, we obtain 

{Id ®Id® Zd) X^od^L {Fd ®Fd®h)M = X^od^div {Fd ®Fd® Id) MZd 

and 

{Id ® /d ® Xd) X^od^L (Fd ®Fd®Id)M = X^od^L (Fd ®Fd® Id) MXd. 
We can rewrite this as 

{Id ®Id®cr)N = Na (19) 
with X = Xmod^div(-^(i ® Fd® Id)M. Hcuce, using the notation of Th. IHlwe have 

A:=X^odZ\i,{Fd^Fd®Id) and B = Id 

since a is an irreducible representation. The matrix X is an element of the intertwin- 
ing space Int(©j^Q^cr, 0"). Following Lemma it has the decomposition 

N = \^i)®Id e C^'^'^ 

with |$i) G <C^^ . Elementary but cumbersome computations'^ show 

\'^i) = {Fl®Id)[Y.\q){q\®X-'n (VdJ2v^,k\j)^\k)] ■ (20) 

We extend the representation a on the right side of Eq. |T9|l to the direct sum of d'^ 
copies of a. The matrix W of the resulting intertwining space has the decomposition 
C ® Id with C G C"'^^'^^. Therefore, we extend {|$i)} to an orthonormal basis 
{|$i), |$2), • • • , I'^'da)} of C'^^ We can define the unitary 

U ■= A^W{B ®B) = At(^(|$i)|$2) . . . \^d^)) ® Ic 



^ Write = X^j^o -^d^j with appropriate diagonal matrices Aj and powers of the shift Xd- 
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that extends M with B := I(n-i)d- Now we show how to simphfy the implementation 
by preparing an appropriate ancilla state. We have 

f/(|0) ® 1^)) = At((|$i)|$2) . . . ® /,)(|0) ® 1^)) = At(|$i) ® 1^)). 

Hence, we can omit the implementation of W if we initialize the ancilla with of 
Eq. ()20|) . In summary, we have to implement the unitary 

(F] ® Fl ® h)Z^^Xl^^{Fl ®h® h) \q) (g| ® X-/ l}j 

after we have initialized the ancillas with the state vector 

d-l 

\l):=VdY,VJ^,,\j)®\k)e€'\ (21) 

j,k=0 

As a special case consider the initial operator fi = \a){a\/d with \a) G and 
{a\a) = 1. In this case we have y/JI = \a){a\/^/d. Furthermore, we have 

I7) = \a) (g) \a). 
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